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Recently, Locatelli and Schoen proposed a transformation of the potential energy that aids the 
global optimization of Lennard- Jones clusters with non-icosahedral global minima. These cases are 
particularly difficult to optimize because the potential energy surface has a double funnel topography 
with the global minimum at the bottom of the narrower funnel. Here we analyse the effect of this type 
of transformation on the topography of the potential energy surface. The transformation, which 
physically corresponds to a compression of the cluster, firstly reduces the number of stationary 
points on the potential energy surface. Secondly, we show that for a 38-atom cluster with a face- 
centred-cubic global minimum the transformation causes the potential energy surface to become 
increasingly dominated by the funnel associated with the global minimum. The transformation has 
been incorporated in the basin-hopping algorithm using a two-phase approach. 



I. INTRODUCTION 

One of the most important types of global optimiza- 
tion problem, and one which is particularly of interest 
to chemical physicists, is the determination of the low- 
est energy configuration of a molecular system, such as a 
protein, a crystal or a cluster .EJ However, such a task can 
be very difficult because of the large number of minima 
that a potential energy surface (PES) can have — it is gen- 
erally expected that the number of jjninima of a system 
will increase exponentially with size.BI Therefore, if appli- 
cations to large systems with realistic descriptions of the 
interatomic interactions are to be feasible, it is necessary 
that efficient global optimization algorithms, which scale 
well with system size, are developed. 

A key part of this development is understanding when 
and why an algorithm is likely to succeed or fail, be- 
cause, as well as providing useful information about the 
limitations of an algorithm, this physical insight might 
be utilised in the design of better algorithms. This is the 
motivation behind the current paper. Here, we analyse 
the reasons for the success of a recent algorithm when 
applied to the global optimization of Lennard- Jones (LJ) 
clusters for some particularly difficult sizes. 

The global optimization of LJ clusters has probably 
become the most comnioB. benchmark for configurational 
optimization problemsw'EI Putative global minima have 
been obtained for all sizes up to 309 atoms,Q liJ and up- 
to-date databases of these structures are maintained on 
the web.BEa There are two types of difficulty for the LJ 
cluster problem. First, there is the general increase in 
the number of minima with cluster sizeHJu Second, on 
top of this effect there are size-specific effects related to 
the topography of the PESj£3 

For most of the clusters the topographsi-ipf. the PES 
aids global optimization. There is a funnelMEJ from the 
high-energy liquid-like clusters to the low etiergy minima 
with structures based up on the MackayO icosahedra. 
When there is a dominant low-energy icosahedral mini- 
mum at the bottom of the funnel, such as when complete 
Mackay icosahedra can be formed, global optimization is 



particularly easy. 

However, there are some sizes for which the global min- 
imum is not icosahedral. At iV=38 the global muiimum is 
a face-centred-cubic (fee) truncated octahedronETEHl (38A 
in Figure 0), at A^=75-7Z-and 102-iq4,tlie global min- 
ima are based on MarksE£l decahedraHEll (e.g. 75A in 
Figure |l]), and at A^=98 the global minimum is a Leary 
tetrahedronllZl (98 A in Figure |]). For these sizes the PES 
has a fundamentally different character. As well as the 
wide funnel leading down to the low-energy icosahedral 
structures, there is a much naixpwer funnel which leads 
down to the global minimum.E30 Relaxation down the 
PES is much more likely to take the system into the wider 
funnel where it is then trapped. liie time scale for in- 
terfunneJ equilibration is-wery slowc3 because of the large 
energyE2l and free energyE3 barriers between the two fun- 
nels. 



34A 38A 75A 98A 




FIG. 1. The global minima and some low- lying minima of 
LJ34, LJ38, LJ75 and LJ98. 34A, 38C, 75C and 98B are based 
on Mackay icosahedra. 34H and 98A are Leary tetrahedra. 
38A is a face-centred-cubic truncated octahedron, and 75A 
is a Marks decahedron. The letter gives the energetic rank 
of the minimum, i.e. global minima are labelled with an 'A', 
etc. 

As a result these eight clusters are hard to optimize, 
the larger examples being virtually impossible to opti- 
mize by traditional approaches, such as simulated an- 
nealing. However, these cases are solvable by a set 
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of methods in which the_ibasin-hopping' transformation 
is appHed to the PES.E3 This traMformation is used 
by the Mmate Carlo minimizationcj or basin-hopping 
algorithm,ll3 andjiopjlicttly by aU the most successful ge- 
netic algorithmsJl3'E2rE3 The transformation of the PES 
works by changing the thermodynamics of the clusters 
such that the system is now able to pass between the fun- 
nels more easilyO'Lj However, the non-icosahedral global 
minima still take miich longer to find than the icosahe- 
dral global minima,Ej and there is no way of knowing if 
one has waited long enough to rule out the possibility 
of a non-icosahedral global minimum. This is illustrated 
by the Leary tetrahedron at N—98. Despite the fact 
that powerftirL-pptimization techniques had been applied 
to LJ98,Ejc3'E3,the global minimum was discovered only 
very recentlyO Subsequently, it was confirmed that this 
minimum coul d . be . found by some of the previously ap- 
plied methods.caeJ 

Given this background, it would be useful to develop 
techniques that are more efficient for these double-funnel 
examples. Two potential approaches have very recently 
been put forward. First, Hartke has achieved improve- 
ments in the genetic algorithm approach by forcing the 
system to maintain a diversity of structural types in 
the population, thus preventing the population becoming 
concentrated in the icosahedral funnel.Ej Second, Schoen 
and Locatelli noted that the exceptions to the icosahe- 
dral structural motifs are usually more spherical than the 
competing icosahedral structures. This is because the 
exceptions generally occur at sizes where both a particu- 
larly stable form for the alternative morphology is possi- 
ble and the icosahedral structures involve an incomplete 
overlayer. Therefore, Schoen and Locatelli added a term 
to the potential energy favouring compact clusters. Us- 
ing this PES transformation, the non-icosahedral global 
minima at A^=38, 98,102-104 were much more likely ta 
be found by their multi-start minimization algorithm.cil 
An additional transformation had to be applied in order 
to find the global minima at N— 75-77. 

It is the reasons for the success of this second approach 
that we examine in this paper. In particular we show how 
Schoen and Locatelli's transformation affects the topog- 
raphy of the PES. We also show how the transformation 
can be incorporated as an element of an existing algo- 
rithm, namely basin-hopping. 



II. METHODS 



The atoms iri^he clusters interact via the Lennard- 
Jones potential^ 
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term proportional to J2i<j '''ij which penalizes long pair 
distances.0 Here, we use a slightly different form, which 
again acts to compress the cluster. The energy for such 
a compressed Lennard- Jones (CLJ) cluster is given by 



^-CLJ — E 
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where /icomp is a parameter that determines the magni- 
tude of the compression acting on the cluster, and rc.o.m 
is the position of the centre of mass of the cluster. We 
found the additional term to be approximately propor- 
tional to Schoen and Locatelli's expression, and so the 
effect of the two transformations on the PES topography 
are virtually identical. 

To map the PES topography of these CLJ clusteijs-iwe 
use the same methods as those we have applied to LjE3'EII 
and Morsec3 clusters to obtain large samples of connected 
minima and transition states that provide good repre- 
sentations of the low-energy regions of the PES. The 
approach-involves repeated applications of eigenvector- 
foUowingQ to find new transition states and the minima 
they connect. r— in 

In the basin-hopping algorithm,E3L3 the transformed 
potential energy is given by 



£;(x) = min{S(x)} , 



(3) 



where e is the pair well depth and 2^/^ct is the equilibrium 
pair separation. To this, Schoen and Locatelli added a 



where x represents the vector of nuclear coordinates and 
min signifies that an energy minimization is performed 
starting from x. Hence the energy at any point in con- 
figuration space is assigned to that of the local minimum 
obtained by the minimization, and the transformed PES 
consists of a set of plateaus or steps each corresponding 
to the basin of attraction surrounding a minimum on the 
original PES. This PES is then searched by constant tem- 
perature Monte Carlo. Additionally, the algorithm has 
been found to be more efficient for clusters if the config- 
uration is rese|ti;o that of the new local minimum at each 
accepted stepo 

There are two ways that one might incorporate a fur- 
ther PES transformation into this algorithm. One could 
use basin-hopping to first find the global minimum of 
the transformed PES, then reoptimize the niow lowest 
energy minima under the original potential. However, if 
the global minimum of the original PES is not among the 
niow lowest energy minima of the transformed PES this 
approach is bound to fail. 

Alternatively, at each step one could first optimize a 
new configuration using the transformed potential, then 
reoptimize the resulting minimum using the original po- 
tential. By incorporating this second minimization the 
shortcomings of the first approach are avoided. Further- 
more, if the energy of this final minimum is used in the 
Metropolis acceptance criterion, the Boltzmann weight 
of each minimum is unchanged. However, the occupation 
probability of a particular minimum will be proportional 
the area of the basin of attraction of the minimum on the 
transformed rather than the original PES, i.e. 
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oc n,Ai exp(-/3i?i), 



where rii is the number of permutational isomers of i 
and Ai is the total area of the basins of attraction of 
the minima on E which when reoptimized on E lead to 
minimum i. Therefore, if the relative area of the global 
minimum is larger on the transformed PES, optimization 
should be easier using this approach. We refer to this ver- 
sion of the basin-hopping algorithm as two-phase basin 
hopping. This variation is not much more computation- 
ally demanding than standard basin-hopping because the 
starting point for the second minimization is likely to be 
close to a minimum of the untransformed PES. 

There is one further difference from previous imple- 
mentations of the basin-hopping algorithm. Previously, 
we had performed minimization in Equation (H) by 
conjugate gradient O However, we have since found pS, 
limited memory BFGS algorithm that is more efficient.cSI 



III. RESULTS 

In global optimization the aim of transforming the po- 
tential energy surface is to make the global minimum eas- 
ier to locate. Typically, one therefore wants the transfor- 
mation to reduce the number of minima and the barriers 
between them. Furthermore, if the transformation is to 
change the relative energies of the minima, one wants the 
energetic bias towards the global minimum to increase. 

As the number of minima and transition states on 
the CLJ13 PES is small enough that virtually all can 
be found, we can examine whether the compressive term 
has the first of the above effects by examining CLJ13 as a 



(4) we obtain a good representation of the lower energy re- 
gions of the PES. At each value of ficomp we obtained 
a sample of 6000 minima. The effect of ^comp on the 
number of stationary points, which we noted for CLJ13, 
is again evident (Table ^). As /icomp increases, risoarch, 
the number of minima from which we have to perform 
transition state searches in order to generate the 6000 
minima, increases and it becomes more likely that a new 
transition state does not connect to a new minimum, but 
rather to one already in our sample. 



function of /^comp ■ The number of minima and transition 
states clearly decreases as //comp increases (Table |). It is 
interesting to note that minima with low symmetry pref- 
erentially disappear. The PES transformation places the 
cluster in a harmonic potential about its centre of mass. 
This potential plays a role similar to a soft spherical box, 
and so less compact minima disappear from the PES as 
/^comp increases. Similar results are found when periodic 
boundary conditions are applied — the number of minima 
is much less than for a LJ cluster of equivalent size and 
the number of minima decreases as the pressure in the 
cell is increased. Cj'lHI 

It is also worth noting that the magnitude of the down- 
hill barriers relative to the energy difference between the 
minima decreases as fj-comp increases (TableJn). In the 
terminology used by Berry and coworkers,E3 the pro- 
files of the pathways to the global minimum become 
more staircase-like and less sawtooth-like with increasing 
/^comp- The combination of the changes to the number 
of stationary points and the barrier heights act to make 
relaxation to the icosahedral global minimum easier as 
the PES is further transformed. 

Next, we examine the CLJ38 cluster. For a cluster of 
this size it is not feasible to obtain a complete representa- 
tion of the PES in terms of stationary points, so instead 




FIG. 2. Qcomp for the LJjv global minima. To make 
the size-dependence more clear the zero is taken to be the 



function, Qave, a four paramater fit to the Qc 



values. 



comp where 
l^/cr^ is 



Qave = 25.915iV - 166.956iV^/^ + 382.765iV^/^ - 293.972. 
Also included in the figure are isolated data points (crosses) 
corresponding to the non-global minima illustrated in Figure 
and the second lowest energy minima for N=76, 77 and 
102-104. 

The second desired effect of a PES transformation is to 
change the energetics in a manner that makes the global 
minimum more favourable. We can get a simple guide as 
to how the energies of the minima depend on /^comp if we 
assume there is no structural relaxation in response to 
changing /icomp- Then Eclj = -Elj + Mcomp 
the order parameter, Qcomp = J2i ~ ^ 
evaluated at /Xcomp=0. From the values of Qcomp we can 
predict the changes in the relative energies of any two 
minima. 

Qcomp is a measure of the compactness of the cluster, 
and from Figure |^ one can see how the compactness of the 
global minima depends on size. For the first two shells 
the icosahedral global minima are most compact when 
complete Mackay icosahedra can be formed, e.g. iV=13 
and 55. However, for the third shell the most compact 
icosahedral structure is at A'^=135, where twelve vertex 
atoms of the Mackay icosahedron are missing, rather than 
at N=U7. 

If we examine LJ38 as an example of a cluster with 
a non-icosahedral global minimum, we see that this size 
corresponds to a pronounced minimum in Figure |— the 
truncated octahedron is particularly compact compared 
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to the other global minima of similar size. Further- 
more, from Figure ||b we can see that the LJ38 global 
minimum has the lowest value of Qcomp of all the LJag 
minima. Therefore, the energy gap between the global 
minimum and the lowest-en erg y icosahedral minimum 
increases with /icomp (Table |ll|). To visualize how this 
deepening of the fee funnel changes the PES topography 
we present disconnectivity graphs of CLJ38 for a range 
of /Xcomp values in Figure ^. 

Disconnectivity graphs provide a repr esentation of the 
barriers between minima on a PES.EBIEj In a disconnec- 
tivity graph, each line ends at the energy of a minimum. 
At a series of equally-spaced energy levels we compute 
which (sets of) minima are connected by paths that never 
exceed that energy. We then join up the lines in the dis- 
connectivity graph at the energy level where the corre- 
sponding (sets of) minima first become connected. In a 
disconnectivity graph an ideal single-funnel PES would 
be represented by a single dominant stem associated with 



the global minimum to which the other minima directly 
join. For a multiple-funnel PES there would be a number 
of major stems which only join at high energy. 

From the disconnectivity graph of LJ38 one can de- 
duce that the cluster has a double-funnel PES (Figure 
^). There is a narrow funnel associated with the global 
minimum, and a wider funnel associated with the icosa- 
hedral minima. There are a number of low-energy min- 
ima at the bottom of the icosahedral funnel, which, al- 
though they have only small differences in the way the 
outer layer is arranged (e.g. the second lowest icosahedral 
minimum, 38C, is depicted in Figure 0), can be separated 
by moderate-sized barriers. As a result there is a certain 
amount of fine structure at the bottom of the icosahedral 
funnel with not all minima joined directly to the stem of 
the lowest-energy icosahedral minimum. From the data 
in Table || one can see that there are many more minima 
associated with the icosahedral funnel. 
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Energy /e Energy /£ 

FIG. 3. Scatter plots of Qcomp against minimum energy for large samples of minima for (a) LJ34, (a) LJ38, (a) LJ75 and 

(d) LJgg. The minima depicted in Figure |l]are labelled by the letter corresponding to their energetic rank. In (d) there are 

two subsets: diamonds correspond to minima found when the search was started from the tetrahedral global minimum and 

crosses correspond to the set when started from the lowest energy icosahedral minimum. There is no overlap between these 

two sets because no pathways connecting the two funnels were located. The patterns of points for Qunoar ~ X^i<j ^''^ 

virtually identical to those of this figure, showing that the current transformation is effectively equivalent to Locatelli and 

Schoen's. 
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FIG. 4. Disconnectivity graphs of CLJ38 for /icomp = (a) (b) 0.25e (c) l.Oe and (d) 5e. In (a) the 150 lowest-energy 
minima are represented in the graph, and in (b)-(d) the 250 lowest-energy minima are represented. The icosahedral and 
fee funnels are labelled. The units of energy on the vertical axis are e. 
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FIG. 5. Heat capacity curves for CLJag with the values 
of /.tcomp/e, as labeUed. The curves were calculated from our 
samples i-eifj-fiOOO minima using the harmonic superposition 
method. EJ'l3 

As the fee funnel beeomes deeper with increasing /icomp 
it increases in size relative to the icosahedral funnel (Fig- 
ure ||). By ^(,oinp=5e the fee funnel dominates the PES, 
and the disconnectivity graph has the form expected for 
an ideal single funnel with only a very small sub-funnel 
for the icosahedral minima. These changes are also re- 
flected in the number of minima associated with both 
funnels (Table |l|. 

These changes to the PES topography of course affect 
the thermodynamics. For LJ38 there are two peaks in 
the heat capacity curve (Figure H). The first is due to 
a transition from the fee global minimum to the icosa- 
hedral minima, which is driven by the greater entropy 
of the latter. The second corresponds to melting. The 
first transition hinders global optimization because it is 
thermodynamically favourable for the cluster to enter 
the icosahedral funnel on cooUgg-from the molten state, 
where it can then be trapped.EHH However, as /icomp in- 
creases, the decreasing entropy of the icosahedral funnel 
can no longer overcome the increasing energy difference 
between the global minimum and the icosahedral funnel 
(Table H) and so this first transition is suppressed. Con- 
sequently, the heat capacity curves for the CLJ38 clusters 
in Figure |^ show only one peak, indicating that the global 
minimum is most stable up to melting. 

Of course, the changes to the PES topography and 
thermodynamics mean that on relaxation down the PES 
the system is more likely to enter the fee funnel as /icomp 
increases. Furthermore, the energy barrier to escape 
from the icosahedral funnel relative to the energy dif- 
ference between the bottoms of the two funnels becomes 
smaller (Table thus making escape from the icosahe- 
dral funneLeasier. To quantify these effects we performed 
annealingEj simulations for CLJ38 at a number of values 
of Hcomp (Table For LJ38 80% of the longer an- 

nealing runs ended at the bottom of the icosahedral fun- 
nel, and only 2% at the global minimum. However, by 
Mcomp=5e 99.5% of the long annealing runs reached the 



global minimum. 

Given the above, it is unsurprising that two-phase 
basin-hopping finds the global minimum more rapidly as 
Mcomp increases (Figure gb). At large /icomp the first- 
passage time is 40 times shorter than for LJ38. Con- 
versely, the first-passage time to reach the icosahedral 
minimum 38C increases. These changes are driven by 
changes to Ai in Equation The basin of attraction 
of the global minimum increases in size relative to those 
of the icosahedral minima as the PES is further trans- 
formed. 



(a) 




34H 



i^. 



(b), 




FIG. 6. The /icomp-dependence of the first-passage time 
(in MC steps) to find the specified minima of (a) LJ34 and 
(a) LJ38 from a random starting configuration in two-phase 
basin-hopping runs. Each point represents an average over 
400 runs. The temperature used is l.Oefc"^. 

Locatelli and Schoen's transformation works for LJ38 
because the global minimum is the most compact spheri- 
cal minimum. However, this does not necessarily have to 
be the case, even for those clusters with non-icosahedral 
global minima. From Figure ^ one can see that the non- 
icosahedral global minima at 7V=98 and 102-104 have 
particularly low values of Qcomp and Figure ^ confirms 
that the Leary tetrahedron, 98A, has the lowest Qcomp 
value of all the LJ98 minima. Therefore, Locatelli and 
Schoen were able to locate these global minima. How- 
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ever, for N=75-77 the values of Qcomp for the Marks 
decahedra are not set apart from the nearby icosahedral 
global minima (Figure ||) and Figure ^ shows that there 
are a number of LJ75 minima which have lower values of 
Qcomp than 75A. In particular, the icosahedral minimum 
75C that is third lowest in energy has a lower Qcomp, 
and the Marks decahedron is no longer the CLJ75 global 
minimum beyond /icomp=3.1e. 

The geometric root of this behaviour is that the Marks 
decahedra at N =7^-77 are the least spherical of the non- 
icosahedral global minima. The 75-atom Marks deca- 
hedron is somewhat oblate and some of the icosahedral 
minima with which 75A is competing are prolate by a 
similar degree, leading to comparable values of Qcomp- 
Therefore, although the transformation may aid global 
optimization by reducing the number of minima and by 
increasing the energy of many minima relative to the 
Marks decahedron, unlike for LJ38 it does not remove the 
fundamental double- funnel character of the PES. To lo- 
cate the global minimum Locatelli and Schoen had to add 
an additional 'diameter penalization' to the potential.tiil 

Locatelli and Schoen found that for many of the clus- 
ters their transformation did not aid global optimization. 
This was not unexpected, but simply reflects the fact that 
often the icosahedral global minima are not the most 
compact minima. We analyse one example. At 7V=34 it 
is possible to form a compact Leary tetrahedron (34H in 
Figure |l|), which is the eighth lowest-energy LJ34 mini- 
mum. This structure has a significantly lower value of 
Qcomp than the global minimum (Figure ^). As a re- 
sult, the Leary tetrahedron becomes the CLJ34 global 
minimum at /Xcomp=0.3e. The results of two-phase basin- 
hopping runs are similar to those for LJsg in that as ^com-p 
increases the compact non-icosahedral structure becomes 
significantly easier to locate and the low-energy icosahe- 
dral minima more difficult (Figure |6|). The difference, 
though, is that now this scenario is undesirable, because 
it is the global minimum that is becoming more difhcult 
to reach. 



IV. CONCLUSIONS 

By analysing the effect of a compressive transforma- 
tion on the PES topography we have obtained insights 
into the reasons for its success in aiding the optimization 
of LJ clusters that have non-icosahedral global minima. 
Firstly, we have shown that the transformation reduces 
the number of minima and transition states on the PES. 
Secondly, for examples where, as is often the case, the 
non-icosahedral global minimum is the most compact 
structure, the transformation causes the funnel of the 
global minimum to become increasingly dominant. For 
LJ38 the PES has a double funnel, whilst at large /icomp 
the PES has an ideal single-funnel topography, enabling 
the system to relax easily down the PES to the global 
minima. However, when as for LJ75, the decahedral 



global minimum is only one of the more compact min- 
ima, the transformation is less beneficial for global opti- 
mization. By contrast, for sizes with icosahedral global 
minima the transition is often unhelpful, as we saw for 
LJ34, because the global minimum is much less likely to 
be the most compact structure. Therefore, the trans- 
formation needs to be used in combination with other 
methods. As the transformation is most likely to be suc- 
cessful for clusters where other methods fail it can act 
as a good complement to them. For example, when the 
basin-hopping algorithm is applied, usually a series of 
runs are performed at each size. If one of the runs used 
the two-phase approach, this would increase the chance 
of success for those sizes where the PES had a multiple- 
funnel topography. 

Other PES transformations could also be usefully em- 
ployed alongside standard basin-hopping runs in this two- 
phase approach, if they are likely to aid global optimiza- 
tion for some sizes. For example, increasing the range 
of the potential is another transformation tliat reduces 
the number of stationary points on the PESJlj Using the 
transformations alongside standard runs avoids one of 
the major difficulties associated with PES transforma- 
tions. They are rarely universally effective, but rather 
there are likely to be some instances when they destabi- 
lize the global minimum, thus making optimization more 
difficult. This is certainly the case when increasing the 
range of the potential, where the range-dependence of the 
most stable cluster structure is well-documented.E2rEll 

Although we have seen how a compressive transfor- 
mation can be useful in aiding the global optimization 
of LJ clusters, an important question is how generally 
useful it will be. Although this question can only be 
definitively answered through applications to a variety 
of systems, one would expect it to be useful for metal 
and simple molecular clusters that form compact struc- 
tures, particularly those that favour 12-coordination. For 
these systems, as with LJ clusters, the strength of this 
approach would be locating those global minima that are 
not based on the dominant morphology, because the al- 
ternative morphologies are only likely to be most stable 
when they are compact and sherical. It might also be 
useful in systems such as proteins where there are a large 
number of less compact unfolded configurations. How- 
ever, it would not be useful for clusters of substances, 
such as water and silicon, which form open network struc- 
tures where the liquid can be denser than the solid. 
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TABLE I. The number of minima, rimin, and transi- 
tion states, nts, for CLJ13 as a function of /Ucomp- For 
each minimum 30 transition state searches were performed; 
these searches were parallel and antiparallel to the eigen- 
vectors with the fifteen lowest eigenvalues. AS, b„, ba 
are the average energy difference, uphill barrier and down- 
hill barrier, respectively, where the average is over all the 
non-degenerate rearrangement pathways. (Degenerate path- 
ways connect different permutational isomers of the same 
minimum.) AE = bu — ba- 
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AE/e 1.593 3.172 4.501 7.191 11.215 20.701 40.176 

bu/e 2.201 3.939 5.396 8.231 12.346 21.759 42.263 

bd/e 0.609 0.767 0.896 1.041 1.131 1.058 2.087 

bd/AE 0.382 0.242 0.199 0.145 0.101 0.051 0.052 



TABLE IL Properties of the CLJgg PES for a sample of 
6000 connected minima as a function of //comp- nts is the 
number of transition states connecting these minima. A_B is 
the energy difference between the global minimum and the 
lowest energy icosahedral minimum and fefcc (&icos) is the en- 
ergy barrier that has to be overcome to escape from the fee 
(icosahedral) funnel and enter the icosahedral (fee) funnel. 
Of course, AE = 6fcc — hcos- For the nscarch lowest-energy 
minima 20 transition state searches were performed; these 
searches were parallel and antiparallel to the eigenvectors 
with the ten lowest eigenvalues, nfcc and the num- 

bers of minima in the fee and icosahedral funnels at the en- 
ergy at which the two funnels become connected. 
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TABLE IIL Results of annealing simulations for CLJ38 as 
a function of /icomp- ./o,, ("cycles) is the fraction of the anneal- 
ing runs that terminated at the global minimum, and /icos is 
the fraction of runs that ended in the lowest five icosahedral 
minimum. Each annealing run involves a linear decrease in 
the temperature from the liquid to OK in ncycies Monte Carlo 
cycles. The results are averages over 200 annealing runs. 
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